Low-voltage current noise in long quantum SINIS junctions 
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The current noise in long SINIS junctions at low temperatures is sensitive to the population of the 
sub-gap states which is far from equilibrium even at low bias voltages. Nonequilibrium distribution 
establishes due to an interplay between voltage-driven inter-level Landau-Zener transitions and intra- 
level inelastic relaxation. We show that the Fano factor is enhanced drastically, being proportional 
to the number of times which particle can fly along the Andreev trajectory before it escapes from 
the level due to inelastic scattering. Combining the dc current and noise measurements one can 
fully characterize the non-equilibrium kinetics in SINIS junctions. 

PACS numbers: 73.23.-b, 74.45. +c, 74.78.Na 



Charge transport through mesoscopic normal- 
superconductor (NS) hybrid structures is a subject of 
intense current interest (see P, 0] for a review). The 
multiple Andreev reflections (MAR) at the NS interfaces 
lead to several unusual features, in particular to the 
so-called sub-gap structure in the current-voltage (I-V) 
curves 0, 0, [If • The electric current noise in NS hybrid 
systems provides information complementary to the 
dc-studies and exhibits properties that are essentially 
different from those of the contacts between the normal 
metals. In normal contacts, the shot noise induced by 
the Pauli correlation is relatively weak Q and vanishes 
if the interface is ideal. In NS systems, on the contrary, 
Andreev scattering enhances reflection of electrons and 
results in a finite shot noise even when the contact itself 
is noiseless in the normal state Q. The situation at low 
voltages is especially interesting: the interplay between 
the quasiparticle scattering and the superconducting 
coherence reveals an intricate structure of the non- 
equilibrium population of sub-gap states and leads to 
considerable enhancement of the noise [H, 0, QjJ II , HJ ■ 



In the present work we study the noise in a single-mode 
long SINIS junction (I stands for an insulator) where the 
length d of the normal conductor is much larger than the 
superconducting coherence length, but shorter than the 
inelastic mean-free path. This is an exemplary hybrid 
structure [l3| where non-equilibrium dynamics is gov- 
erned by inelastic scattering and voltage-driven Landau- 
Zener (LZ) transitions between Andreev sub-gap states. 
In particular, in a long junction where the number of the 
states involved in the transport is large, the LZ transi- 
tions lead to a significant increase in the dc current far 
above the critical Josephson current of the junction 14\. 



and is related to thermal fluctuations in the occupancies 
of the sub-gap states [1, 0, [l(| ■ The width of the zero- 
frequency peak in the noise spectrum is determined by 
the inelastic relaxation rate and thus provides an exper- 
imental tool for measuring its temperature dependence. 
We show that for low temperatures and small bias volt- 
ages the noise is controlled by the non-equilibrium popu- 
lation of the sub-gap states characterized by an effective 
temperature T e ff which is determined by the competition 
between the inelastic scattering and the voltage-driven 
inter-level LZ transitions; T e ff depends on the bias volt- 
age and remains finite even at zero ambient temperature. 
At the same time, the dc current depends not only on T e g 
but also on a LZ-controlled parameter that accounts for 
the charge transfer from the excited levels. For weak 
LZ processes, this parameter and thus the dc current are 
small. As a result, the ratio of the zero-frequency noise 
to the dc current, i. e., the Fano factor proportional to 
the fluctuating "effective charge" is increased drastically 
by a factor qualitatively different from that predicted in 
Refs. 0, [l(| for short contacts. Since both the dc I-V 
curve and the noise spectrum depend on inelastic relax- 
ation rate and on the LZ probabilities, the combined I-V 
and noise measurements enable detailed characterization 
of non-equilibrium dynamics in long SINIS junctions. 
Noise spectrum. The current operator is 



f. ieh 

I(ti,xi) = 7T- 
2m 



(Vx-Vi) ^t(i)VUi') 



(1) 



J 1=1' 



The noise properties of such systems are very informa- 
tive. At finite temperature, the current noise is finite 
even in the absence of a dc bias voltage (and dissipation) 



where ip a is the operator of particle field with a spin a. 
Noise spectrum is defined as a Fourier transform, 

S(u) = Je^dr(^Ai(t-^x),Ai(t + ^,x)}) (2) 

where AI = I — (I), (I) is the statistical average current, 
and {A, B} = (AB+BA)/2, while overline means the av- 



2 



erage over a long time interval, To: A = T _1 f^ A(t) dt. 

We focus on the low-frequency behavior of the noise 
correlator, such that Huj is much smaller than the char- 
acteristic energies of the states involved. In this situa- 
tion, making use of the Bogoliubov transformation one 
can express the correlator in Eq. ([2]) through the occu- 
pation numbers /(e„) of the sub-gap states and matrix 
elements, L mn , quadratic in current operators: 



4m 2 



■Jj5(e«m + hhj)L nm [l - /i(e„)/i(e m )] , (3) 



L nm (x,t) = {(Vi - Vi) (V 2 - V 2 ) [<(l)n TO (l') (4) 
+v* n (l)v m (l')] K(2')<(2) +i; n (2')<(2)]} 1=1 , =2=2 , . 

Here u n , v n and e„ are the wave functions and energies 
of the sub-gap states, which have to be found from the 
Bogoliubov-dc Genncs (BdG) equations, e nm — e„ — e m . 
These quantities depend adiabatically on time t. The 
matrix elements L nm are symmetric in the transposition 
n <-> m. The functions f(e n ) are the occupation numbers 
of the subgap states; here we have introduced the odd- 
in-e combination /i(e) = /(— e) — /(e). We will see later 
that the even combination does not appear in our case, 
f 2 = 1 - /(e) - /(-e) = 0. 

Sub-gap states. To find the matrix elements L mn we 
solve the BdG equations for e < |A| modelling the SI- 
NIS structure as two insulating barriers described by a 
symmetric 5-function potential with a magnitude /. By 
adjusting the barrier height / one can account both for 
the velocity mismatch in the adjacent conductors and for 
the tunneling of electrons through the interface. The gap 
A = \A\e iXR - L in the right, x > d/2, and left, x < -d/2, 
superconductors has xb„l = respectively, while 

A = in the normal conductor, —d/2 < x < d/2. 
The states depend on the phase parameters a = k x d 
and (3 = ed/hv Xl as well as on the unitary scattering 
matrices [l| 5 at each NIS interface 

S N e lS S A e ix 



5 = 



Sa& lx Sne 



-id 



(5) 



which is parameterized through the barrier strength Z = 
ml/h 2 k x and the scattering phase S = arctanl/Z. Here 
k x and v x are the electron wave vector and velocity, re- 
spectively. See Refs. 0, HI for more details and expres- 
sions for Sjy, Sa- 

The wave function inside the normal conductor is a 
linear combination of an incident particle and an incident 
hole, each experiencing Andreev and normal reflections 
at the interfaces. For reflections from the right barrier at 
x = d/2 we have 



u> + 5 A e^ /2 v < + 5j V e l5 w < 
j> + Sa^ 12 ^ + S N e~ lS i 



; (6) 



&><<(x) = e**>+* (lY €<<>{x) = e**>-« ( J ) , 



where q± = k x ± e/Hv x . Comparing Eq. ([6]) with the 
similar equation for the left barrier at x = —d/2 we find 



3 -2i/3 



e^5 2 - e 



-lioC q2 
N 



-2bS A S N cos(a' + 4>/2) 



be 



-2i/3 



n 2ia' Q 2 



/2 C 2 



N 



-2aS A S N cos(a' + </>/2) 







(7) 



(8) 



Here a' = a + 5. The wave function Eq. © is normal- 
ized /(M 2 + \v\ 2 )dx = 1. Defining S N = e^\S N \, the 
solvability of Eqs. CO)-© requires [4[ 



\S N \ 2 sin 2 a' + \S A \ 2 cos 2 (</>/2) = sin 2 (/3 + 7 ) 



(9) 



A long contact has d S> hv x /A and a large number of 
levels N ~ dA/hv x 3> 1. The levels are separated by 
minigaps at <f> = 2nk and <f> = 7r(2fc + 1) which disappear 
when a' — tt/2 + nm and a' = 7tto, respectively. 

Equilibrium noise. For frequencies hw much lower 
than the interlevel spacing e nj „+i ~ hv x /d, only diagonal 
elements contribute to Eq. ([3]) due to the S functions, 



-Wkt 



K\ 2 )\S A \' 



ic/>/2 



a nb* n S* N S ' AC 



-i8+iip/2 



Using Eqs. ©-(O for the cases of nearly ballistic and 
weakly transparent contacts we obtain, respectively, 




sin 2 a' 



Putting /x(e) - /< 0) (e) + /i(e) where f^(e) 
tanh(e/2T) is the equilibrium distribution we find 

S(lu) = - 



z^o, \S N \ 
Z -> oo, | 5a | 



0; 
0. 

(o). 



(10) 
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vi Q) h ■ (ii) 



The first two terms in the square brackets in Eq. 
give the equilibrium thermodynamic noise with 

1 - [/i 0) ( £ )] 2 = cosh~ 2 (e/2T). The 5-function in Eq. ^ 
is broadened by inelastic relaxation processes described 
by the rate r _1 corresponding to energies e <C A. Since 
the interlevel spacing e„.„+i 3> H/t while ftw <C £ n ,n+i 
the 6(cu) should be replaced with 

1 



AM = 7T- 



2tt lu 2 + l/4r 2 ' 
For high transparency we get from Eqs. (fT0| - (fTT|) 



5 t h M 



^A(w)^cosh- 



d 2 



n>0 



2T 



(12) 



If hv x /d <C T <c A, the energies e <C A. Replacing the 
sum with the integral (2d/irhv x ) J de we find 



5 t hH = (4,Tv x e 2 /hd)A(iu) 



(13) 



3 



For a low-transparency junction, Z — > oo, one has 
| | 2 = 1/4Z 4 for e <C A, and the noise acquires an addi- 
tional factor (32Z 8 sin 2 2a') -1 as compared to Eq. (fTB")) . 
The current noise in Eq. Q13p is associated with thermal 
fluctuations in the occupancies of the sub-gap states in 
the same manner as obtained in Ref. [|| for short junc- 
tions. However, the spectral density in long junctions is 
linear in temperature which accounts for the number of 
excited levels N e g = 2Td/irhv x contributing to the noise. 

Non-equilibrium distribution in high-transparency 
junctions. A dc current at a finite bias voltage appears 
only due to nonequilibrium quasiparticle distribution. 
Indeed, in equilibrium the current would be just the 
supcrcurrent adiabatically driven by the time-dependent 
phase difference 4> = uijt + 4>o and having zero time- 
average (wj = 2eV/H is the Josephson frequency). In 
a highly-transparent barrier, particles with |e| > |A| 
remain in equilibrium with the heat bath for eV <C A 
because they have sufficient time to escape from the 
double-barrier region and relax in the continuum. How- 
ever, the time depending phase induces inter-level LZ 
transitions across the minigaps near cf> — nk producing 
deviation from equilibrium for the sub-gap states. As 
the phase changes with time, the outermost Andreev 
levels split off from the continuum and capture particles 
at energies e = ±|A|. These quasiparticles can be then 
excited to the other levels by LZ transitions. The inter- 
play between the LZ transitions and the weak inelastic 
relaxation produces the nonequilibrium population of 
the subgap states and leads to a finite dc current. 

The kinetic equation for the subgap distribution was 
recently derived in Ref. 14| under assumption that LZ 
transitions occur only near the avoided crossings at <j) = 
•nk. This can be realized if the minigaps as well the bias, 
eV , are much less than the typical distance Hv x /d be- 
tween the levels. Only in this case the avoided crossing 
points are well defined and there are no voltage-induced 
inter-level transitions elsewhere. These requirements can 
be met only in almost ballistic contacts, 1 — T <C 1, where 
T is the transparency of the interface. The kinetic equa- 
tion also assumes incoherent LZ transitions. Indeed, the 
inequality eV <C h~v x /d ensures that the correlation time 
of the quantum beating between the adjacent Andreev 
states, ft/|e„ + i. n | w d/v Xl is much shorter than the time 



H/2eV between successive level crossings. 



Between the LZ tunnelling events, the distribution 
function relaxes due to inelastic scattering with the rate 
t , which is assumed to be energy-independent. Since 
the even component is absent, fa = 0, we have /i(e„) = 
1 — 2/(e„). For low temperatures, T <C hv x /d, the equi- 
librium distribution /■£ (e n ) = sign(e„). As a result, 



/iK(0)]-/i (O) (en) 



p<ir, 
< 2n. 



The LZ transitions at (f) = 27rfc and <p = 71 + 27r/c 
with (constant) probabilities po andp^, respectively, lead 



to difference equations for ip n and Xn with constant 
coefficients [LJ. The boundary conditions require that, 
for the levels moving upwards with the increase in <j), the 
distribution coincides with the equilibrium at e = — |A|, 
while, for levels moving downwards, the distribution co- 
incides with the equilibrium at e = +|A|. Solution to 
the difference equations [l4[ obeys ip-n = —ip n +i, which 
implies the property /i(e) = — e) used earlier. For 
e n > we have ip n+1 +ip n = —2e"F(N—ri)/F(N), where 
v = 7rh/2eVr is the degree of inelastic relaxation during 
the Josephson half-period, 

F(z) = e rz (l + e r w + )-e~ rz (l + e- r w_) , 

2poPtt cosh(r) 



w± 



p e ±r + p^e Tr 



C+P0Ptt(1 



3 ±2r 



) + PO + Prr - 2pQPn 



and C, = e? v — 1. Depending on voltage, v can be either 
large or small. For very low voltage, v 3> 1, the relax- 
ation rate is r ~ v, and the distribution relaxes quickly. 
The most interesting case is when the degree of inelastic 
relaxation is small, In this limit, r is given by 



sinh 2 (r) = v{v + p + p„ - 2p p„)/p p w 



(14) 



The inverse rate r 1 describes the energy broadening of 
the distribution and determines an effective temperature 



eff 



r 1 (de n /dn) = {nhv x /2rd) 



(15) 



which can much exceed the inter level spacing if r <C 1. 

Let N = 2dA/irhv x be the number of sub-gap levels 
with e„ > 0. For cases (z) 1/N <C r <C 1 and (ii) r <^ 
l/N we have 



4'n+l + 4'n = -2 < N +(p Q +p 7T - 2p p 7r )- 



N + (po +p v - 2p p^) 1 



(ii). 



(16) 

Nonequilibrium noise. The equilibrium noise vanishes 
if T <C hv x /d. For high transparency Z ~ > and slow 
relaxation 1/ < 1 the nonequilibrium part of Eq. pip is 



o 2 2 N / 2 
fc=0 



ln=2fc+l 



Using Eq. (fTB"|) one finds 



S do (uj) = (2Tre<vi:N e s/d')A(uj); 



(17) 

N cS = I N 2[N{p +p v -2p p 7r )}- 1 + l 



2 [N(po+Pi r -2p Q pn)}~ 1 + 1 



(ii). 



N e g is the number of excited levels. In the case (i) the 
noise S n c{u) has the form of Eq. (fT3|) where the real 
temperature is replaced with the effective temperature 
T cff = nhv x N cS /2d defined by Eq. JT3J> . This result 
agrees qualitatively with the low-voltage noise in short 
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contacts obtained in Rcf. [gj if the real temperature there 
is replaced with T e g. The zero-frequency noise is 

S nc (0) = (4e 2 v*/d 2 )N eS T = (8T cS e 2 v x T/nhd) . (18) 

This noise does not vanish at zero real temperature; on 
the contrary, it diverges for r — > oo since N c ff ~ N and 
T c ff ~ A when v — > 0. The effective temperature depends 
on voltage; it decreases exponentially for low voltages and 
grows linearly for higher voltages (see below). 
The dc current calculated in Rcf. [1J| is 



The factor A = v/r = N e gv is found from Eq. {14 



(19) 



A = 



PoP* 



1 + (Po +Pit - 2p p w )/u 



< 1 . 



It is a LZ-controllcd efficiency of the non-equilibrium dis- 
tribution in contributing to the charge transfer through 
the junction; A = 1 if po = p-n = 1, while A = if 
Po — Pit = 0: there is no dc current if A — > 0. 

One can write 5 nc (0) = 2ql = 2eIF where F is the 
Fano factor, and q = eF is the effective charge q = 2en/X. 
Here n = v x r/d is the number of times which a particle 
can fly back and forth along the Andreev level trajectory 
before it escapes due to inelastic scattering. The effec- 
tive charge is huge g>e since r ^> d/v x . The enhance- 
ment of fluctuating charge for MAR in short junctions 
has been predicted theoretically in Refs. [1, E(| and con- 
firmed experimentally in Refs. [HI, E3 ■ For a finite bias 
voltage the charge is enhanced by a factor n = A/eV 
which is the number of flights of a particle in the nor- 
mal region before it accumulates the energy sufficient to 
escape into continuum. Here we predict that, in long 
junctions, the fluctuating charge is increased by a larger 
factor n/X = v x r/dX for bias voltages eV <C Hv x /d. 

For highly-transparent contacts 1 - T « 1 and low 
energies, the probabilities po and of LZ tunneling 
can be easily calculated. Assuming that (f> = ujjt + cf>Q 
one finds [14j p^ = exp [— (uiq/luj) sin 2 a'] where luq = 
irv x (l — T 2 )/Td. For po one replaces sin a' with cos a'. 
Therefore, for uij <C the efficiency factor A is small, 
leading to enormous enhancement of the effective charge. 
For loj ^> u)q one obtains po + p v — 2pop„ = loq/ujj. In 
the case Nr 3> 1 the effective temperature is linear in 
voltage, T cff = (Hv x /2d)(ujjT/ sji + cj t/tt). The LZ effi- 
ciency A = 1 / yl + lo^t/tt is much smaller than unity for 
u>ot 3> 1. As a result, the effective charge is still much 
larger than 2en, namely q = 2enyT + ujqt/tt. 

In conclusion, we have considered both thermal and 
non-equilibrium low-voltage noise spectrum in long SI- 
NIS junctions. The spectrum has a pronounced peak at 



uj = 0; its intensity diverges at zero temperature, while 
the width is equal to the inelastic relaxation rate. The 
spectrum allows one to determine the effective tempera- 
ture of the sub-gap excitations which remains finite even 
at zero ambient temperature. At weak inelastic scat- 
tering the Fano factor is large; it is determined by LZ 
tunnelling between the sub-gap states. 
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